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Abstract 

According to a first-principles approach, we clarify electric cloaks' universality concerning a new class 
of phase transitions between negative pathway (NGP) and normal pathway of electric displacement fields, 
which are driven by the geometric shape of the cloak. We report that the NGP arises from shape-enhanced 
strong negative electric polarization, and that it is related to a symmetric oscillation of the paired electric 
permittivities, which are shown to satisfy a sum rule. The NGP does not occur for a spherical cloak, but 
appears up to maximum as the ratio a/6 between the long- and short- principal axis of the spheroidal cloak is 
about 5/2, and eventually disappears as a/6 becomes large enough corresponding to a rod-like shape. Then, 
the cloaking efficiency is compared between different geometrical shapes. The possibility of experiments 
is discussed. This work has relevance to crucial control of electric fields and to general physics of phase 
transitions. 

PACS numbers: 41.20.-q, 05.70.Fh 



* Electronic address: 052019029@fudan.edu.cn 

^ Corresponding author. Electronic address: jphuang@fudan.edu.cn 

Typeset by REVTeX 1 



1 Introduction 

Metamaterials 120 are a new class of electromagnetic materials, which are currently under 
extensive investigation Jj, . They owe their properties to subwavelength details 

of structure rather than to their chemical composition, and can be designed to have properties 
difficult or impossible to find in nature. An electronaagnetic spherical cloak is well-known for 
hiding objects from electro-magnetic waves js, I4, 5, 



lOh . due to the high freedom of design of 



metamaterials [l6D, which maybe both inhomogeneous and anisotropic in their electric permittivity 



and magnetic permeability 



Owing to intriguing potential applications, it 
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has received extensive attention, e.g., ranging from its two-dimensional counterpart [|4J, |7|, |80 and 
scattering cross section to its extensions like acoustic cloaks [H]. In this work, by using a 
first-principles approach based on the coordinate transformation method, we investigate electric 
cloaks and clarify their universality concerning a new class of phase transition between negative 
pathway (NGP) and normal pathway (NMP) of electric displacement fields as geometric shape of 
the cloak changes in a typical range. The shape-driven phase transition seems to be of particular 
interest when compared with abundant phase transitions already known in condensed matter and 
nature. We shall also calculate cloaking efficiency of the electric cloaks with various shapes. This 
work has relevance to crucial control of electric fields and to general physics of phase transitions. 

The remainder of this work is organized as follows. In Section 2, based on the coordinate 
transformation method, we present a first-principles approach to derive the expressions for electric 
permittivities in electric cloaks with non-spherical shapes. This is followed by Section 3, in which 
we numerically clarify their universality concerning the phase transition between negative pathway 
and normal pathway of electric displacement fields, and calculate cloaking efficiency under various 
conditions. This paper ends with a discussion and conclusion in Section 4. 

2 Formalism 

To proceed, we shall use a first-principles approach, namely, the coordinate transformation 
method [|3|,|4t]. We consider that each point in the original Cartesian coordinates can be represented 
as (x, y, z), while in the distorted coordinates it can be given as {u, v, w), which is the location 
of the new point with respect to the x, y, and z axes, namely, u = u{x, y, z), v = v{x, y, z), and 
w = w{x, y, z) [3]. In the distorted coordinates, the form of Maxwell's Equations are invariant, but 



the corresponding electric permittivity are got in the following form e'^ = 



(e„: electric 
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permittivity in the original Cartesian coordinates), where 



dv J \dv J \dv ^ 

Qi = {I^X^(^X^(I^X- (3) 



dw J \dw J \dw J 
Now we consider an eUipsoidal cloak in three dimensions. In the Cartesian coordinates (x, y, z), 
the equation for describing an ellipsoidal shape is 

h — H = 1, 

6^ 

where a, 6, and c are the three principal semi-axes of the ellipsoid, respectively. Next, we choose 
(A, /i, u) to represent a point in the ellipsoidal coordinates. So far, the relation between (A, /i, i/) 
and (x, z) is given as 

2 _ (g' + A) (g^ + /i) (g^ + I/) 

(g2 - 62) (^2 - c2) ' ^ ^ 

^2 ^ + A) (6^ + /.) (6^ + z.) 

^ (62 - c2) (62 _ ^2) ' 

2 ^ (c^ + A) (c^ + (c^ + ^) ... 

(c2 - a2) (c2 - 62) • ^ ^ 

Then, we squeeze the ellipsoidal volume into an ellipsoidal shell through the following relations 

\' . ^g2-gi 

A = gi + A , (7) 

fi'^h + i,^-^, (8) 

u' = c, + u'-^, (9) 

C2 

where gi, 6i, and ci (g2, 62, and C2) are the inner (outer) three principal semi-axes of the ellip- 
soidal shell. In the distorted ellipsoidal coordinates, we have the renormalized values of electric 
permittivities, e^, e'^,, and e'^,, inside the ellisoidal cloak (shell): 

_ g2 h2 C2 {\'-aif ^ 

— 1 r 77^ ^host's \^^) 

a2 - gi 62 - 61 C2 - ci 
_ a2 62 C2 {fi'-bi)\ 

^n' ~ 7 T ^host^ \^^) 

^ g2 - gi 62 - 61 C2 - ci 11'^ 
_ a2 62 C2 {i^'-cif 
a2 - gi 62 - 61 C2 - ci z/'2 
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where Chost denotes the electric permittivity outside the cloak. From equations ([T0l)-([T2l). we can 
see that the electric permittivities are both materially anisotropic and spatially inhomegeneous in 
the distorted ellipsoidal coordinates. For achieving them, metamaterials can help [i6J. 
3 Numerical results 

3.1 Phase transition between negative pathway and normal pathway of electric displacement 
fields 

In our numerical simulations, there are three domains in the whole system, namely, the inner 
domain (domain I), the cloaking or shell domain (domain II), and the outer domain (domain III). 
For convenience, we set the electric permittivity of domain I to have the same value as that Chost 
of domain III, and further take Chost to be equal to the electric permittivity of free space eo. This 
will not affect our results at all because the main interest of this work is on geometrical control of 
electric fields within domain II. The electric permittivity of domain II is determined according to 
equations ([T0l)-([T2l). For all the spherical or non-spherical cloaks discussed in Figs. [HISl we keep 
the volume ratio between domain II and domain I at 7: 1, and the applied electric potentials are set 
to be +3 V and -3 V at the two opposite planes of the cubic box, respectively. It is worth mentioning 
that our results are independent of the actual value of the potentials. As the spherical cloak is 
centro-symmetric, it makes no difference when the electric field is applied in various directions. 
However, for an ellipsoidal cloak, we have to discriminate the direction of applied electric fields 
due to the existence of geometric anisotropy. In our simulations, we apply the electric potentials in 
two opposite planes along the three principal axes, a, 6, and c, of the ellipsoidal cloak, respectively. 
We focus on two types of rotational ellipsoid (namely, spheroid) with a < b = c (oblate spheroid) 
and a > b = c (prolate spheroid). Throughout this work, the set of a, 6, and c will be used to 
denote both three inner principal semi-axes (ai, 6i, and ci) and three outer principal semi-axes 
(a2, 62, and C2) of the ellipsoidal cloak if there are no special instructions. We should also remark 
that our simulation results are independent of the length scale of the electric cloaks. 

The streamlines of Fig. [H represent the pathway of the electric displacement field in the electric 
spherical cloak for the parameters as indicated in the caption. Clearly, the displacement field goes 
around the inner domain, but is limited within the cloak (shell), and eventually returns to its orig- 
inal pathway. In this process, an object inside the inner domain has no effect on the displacement 
field at all. In other words, the object is protected from the invasion of electric fields by using the 
electric cloak. (Its cloaking efficiency will further be investigated in Section 3.2.) 

For a spheroidal cloak, the pathway of the electric displacement field is also illustrated in Fig. [21 
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(oblate spheroid) and Fig. [3] (prolate spheroid). For the oblate spheroidal cloak a < 6 = c, the 
pathway is similar to the spherical cloak discussed in Fig. [T] when the applied electric field is 
directed along the short principal semi-axis a, see Fig.[2l However, the situation becomes different 
for the case of prolate spheroid with a > 6 = c, as shown in Fig. [3] in which the external electric 
field is along the long principal semi-axis a. Figure [3] shows a NGP (negative pathway) of the 
electric displacement field inside the spheroidal cloak. Here, the so-called NGP means that the 
displacement field goes backwards, in contrast to NMP (normal pathway) for which the field goes 
forwards. (All the NGP streamlines shown in Fig. [3] are seemingly limited to enter the inner 
domain. However, this is actually not true. Depending on different streamlines/pathways and/or 
shape parameters, such NGP streamlines can also appear within the cloak domain only.) Owing to 
the symmetry of the prolate spheroid, the NGP zone is strictly symmetrical, which locates close 
to the two opposite sides of the prolate cloak. On the other hand, we have also investigated the 
case of a general ellipsoid with a 7^ 6 7^ c for which the external electric field is directed along 
the longest principal semi-axis, the behavior is generally similar to Fig. [3] (no figures shown here). 
Nevertheless, in this case, the NGP zone is no longer strictly symmetrical, the degree of which 
depends on the actual values of a, 6, and c. 

Now let us dig out the physical origin of the NGP mentioned above. It is known that an electric 
displacement field D is given by D = eoE + P, where E is the external electric field, and P din 
electric polarization. Corresponding to the NGP, D has a negative sign, which is opposite to E 
(Strictly, there is an obtuse angle between D and E in the present case, due to the presence of 
anisotropic and inhomogeneous electric permittivities inside the cloak). Thus, this requires that P 
should be sufficiently negative, satisfying P < — cqE. So far, we can conclude that the appearance 
of the NGP zone results from the shape-enhanced strong negative electric polarization. 

Figure Ua) shows a phase diagram for spheroidal cloak with three principal semi-axes a ^ b = 
c, in the presence of external electric fields along the direction of the principal semi-axis a. For 
this phase diagram, the potential difference Uab (see the caption for its definition) is taken to be 
the order parameter, and the ratio between a and b (a : b) the tuning parameter ranging from to 
3.6. In the typical region of a : 6, a phase transition phenomenon comes to appear. In detail, there 
is an existence of two distinct phases, NMP phase and NGP phase, and a transition region between 
them, in which the NGP is enhanced significantly. Such a phase transition is driven by geometrical 
shape of the cloak. To understand more clearly, in Fig. ffib) we investigate [/^^ as a function of 
a : b ranging from to 20. It is evident that when a : 6 is smaller than or equal to 1 (oblate 
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spheroid or sphere), the electric potential Uab = 0, namely, there is no NGP. As a : 6 > 1 (prolate 
spheroid), the NGP comes to appear, and reach maximum at about a : 6 = 5 : 2. Eventually, it 
tends to disappear as a : 6 is large enough corresponding to a rod-like shape. 

Figure [5] illustrates the paired electric permittivities, 6a and e^, in the NGP streamline within 
the cross section of the spheroidal cloak with a > 6 = c as a function of the distance r (see the 
caption), for a : 6 = (a) 2 : 1, (b) 3 : 1, and (c) 4 : 1. All the points displayed in (a)-(c) have been 
taken in sequence from the starting point to the ending point of the NGP streamline, as indicated by 
the ordered numbers. As a : 6 is given, all the paired sx and (corresponding to the same number) 
appear in a symmetric oscillation trajectory, and satisfy a sum rule: (a) £x + = 2.015 ± 0.003, 
(b) 2.218 ± 0.007, and (c) 2.461 ± 0.006. It is worth noting that the electric permittivity keeps 
unchanged for the cross section (in two dimensions) of our interest. 

3.2 Cloaking efficiency 

We define the cloaking efficiency 77 as the total electric energy of domains II and III over the 
whole system of domains I, II, and III. Apparently, if r] is equal to 100%, the cloaking is perfect. 
In this case, any objects within domain I can be perfectly shielded from external electric fields. In 
other words, the degree of imperfectness of the cloaking is denoted by how r] deviates from 100%. 
For consistence, the applied electric potentials and the electric permittivities of domains I, II, and 
III are set to be the same as those already presented in Section 3.1. In the simulations, the volume 
ratio between domain II and domain I is also kept at 7:1. 

The calculated cloaking efficiencies of the electric cloaks under different conditions are shown 
in Tables I-III. As we can see, the cloaking efficiency of the spherical cloak is 99.58%, which is 
very close to, but not equal to 100% (perfect cloaking). The deviation from 100% arises from 
the small amount of central streamlines that unavoidably pass through domain I. Further, we find 
that the cloaking efficiencies of non-spherical cloaks are always smaller than that of the spherical 
cloak. In detail, for prolate spheroidal (a > 6 = c), oblate spheroidal (a < 6 = c), and ellipsoidal 
(a ^ b ^ c) cloaks, their cloaking efficiencies corresponding to three different principal semi-axes 
have only slight difference and are smaller than that of the spherical cloak. For the three kinds 
of non- spherical cloaks, we can conclude that the cloaking efficiency corresponding to the field 
directed along the longest principal semi-axes is the smallest among those along the three principal 
semi-axes, and that the cloaking efficiency of ellipsoidal (a b ^ c) cloaks is the largest whereas 
that of oblate spheroidal (a < b = c) cloaks is the smallest. 

4 Discussion and conclusion 
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Besides DC electric fields adopted in the above simulations, our results also hold for AC fields 
due to the symmetry of the system. Since they have been obtained by using numerical simulations 
according to a first-principles approach, namely, the coordinate transformation method, they are 
exact without any approximation (despite the error which might occur when we extracted the data, 
e.g., as implied by the sum rules obtained from Fig.O. Regarding an experimental demonstration, 
it would be quite convenient to use metamaterials to realize the electric cloaks with a shape-driven 
phase transition between the NOP and NMP, due to the existing significant achievements in the 
field [5, ^. Further, the fact that our results do not depend on specific length-scales makes the 
experiment more tractable. 

To sum up, we have accurately revealed electric cloaks' universality concerning a new class 
of shape-driven phase transitions between the NOP and NMP of electric displacement fields. We 
have showed that the NGP results from shape-enhanced strong negative electric polarization, and 
that it corresponds to a symmetric oscillation of the paired electric permittivities, which have been 
demonstrated to satisfy a sum rule. The cloaking efficiency has also been calculated for various 
geometrical shapes. 

We thank W. J. Tian for fruitful discussion. This work was supported by the Pujiang Talent 
Project (No. 06PJ14006) of the Shanghai Science and Technology Committee, by the Shanghai 
Education Committee and the Shanghai Education Development Foundation (Shuguang project), 
by the National Natural Science Foundation of China under Grant No. 10604014, and by Chinese 
National Key Basic Research Special Fund under Grant No. 2006CB921706. 
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Figure Captions 

Fig. 1. (color online) Spherical cloak (left panel) and its cross section (right panel). The stream- 
lines denote the pathway of the electric displacement field when the applied electric potentials are 
set to be +3 V and -3 V at the top and bottom plane of the cubic box, respectively. Parameters: 
ri = 0.1 m (inner radius) and r2 = 0.2 m (outer radius). 

Fig. 2. (color online) Oblate spheroidal cloak with three principal semi-axes a, 6, and c sat- 
isfying a < b = c (left panel) and its cross section (right panel). The streamlines denote the 
pathway of the electric displacement field in the direction of the short axis a when the applied 
electric potentials are set to be +3 V and -3 V at the top and bottom plane of the cubic box. Pa- 
rameters: ai = 0.1m and 6i = ci = 0.2 m (three inner principal semi-axes); a2 = 0.2 m and 
62 = C2 = 0.4 m (three outer principal semi-axes). 

Fig. 3. (color online) Same as Fig.[2l but for prolate spheroidal cloak a > b = c and the electric 
displacement field in the direction of the long axis a. Parameters: ai = 0.2 m and 61 = ci = 0.1 m; 
a2 = 0.4 m and 62 = C2 = 0.2 m. 

Fig. 4. (color online) (a) Phase diagram for spheroidal cloak with three principal semi-axes 
a 7^ 6 = c, which represents the existence of two distinct phases, NMP phase and NGP phase, 
and a transition region between them. The potential difference Uab between the starting point A 
and ending point B of the NGP streamline [as also indicated in the inset of (b)] is taken to be the 
order parameter, and the ratio between a and b (a : b) the tuning parameter ranging from to 3.6. 
The three curves correspond to the three incident electric fields, which respectively have a vertical 
distance of 0.86, 0.56, and 0.36 with respect to the principal semi-axis a, see also the inset of (b). 
(b) Uab as a function of a : 6 ranging from to 20, and others are the same as (a). 

Fig. 5. (color online) The paired electric permittivities, ex and e^, in the NGP streamline within 
the cross section of the spheroidal cloak with a > 6 = c as a function of the distance r between 
any point in the NGP streamline and the center of the cloak, for (a) a : 6 = 2 : 1, (b) a : 6 = 3 : 1, 
and (c) a : 6 = 4 : 1. All the 20 points have been taken in sequence from the starting point to the 
ending point of the NGP streamline, as clearly indicated in the corresponding ordered numbers, 
see also the three insets. 
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FIG. 1: 
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FIG. 2 
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FIG. 3: 
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FIG. 5: 
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TABLE I: Cloaking efficiency of spherical and prolate spheroidal (a > b = c) cloaks. The external DC 
electric field is applied along the a, b, and c principal semi-axes, respectively, ri (or r2) denotes the inner 
(or outer) radius of a spherical cloak. 



shape 


direction 


inner 






outer 






V 


sphere 




ri = 0.1m 






r2 = 0.2 m 






99.58% 


spheroid 


a 


ai 0.2 m, bi 


= 0.1m, ci 


= 0.1m 


a2 0.4 m, 62 


= 0.2 m, C2 


= 0.2m 


99.14% 


spheroid 


6 or c 


ai = 0.2 m, bi 


= 0.1m, ci 


= 0.1m 


a2 = 0.4 m, 62 


= 0.2 m, C2 


= 0.2 m 


99.15% 
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TABLE II: Cloaking efficiency of oblate spheroidal (a < b = c) cloaks. 



shape 


direction 


inner 








outer 




rj 


spheroid 


a 


ai = 


0.1m, bi 


= 0.2 m, ci 


= 0.2 m 


a2 = 0.2 m, &2 


= 0.4 m, C2 


= 0.4m 97.75% 


spheroid 


& or c 


ai = 


0.1m, bi 


= 0.2 m, ci 


= 0.2 m 


a2 = 0.2 m, 62 


= 0.4 m, C2 


= 0.4 m 97.69% 
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TABLE III: Cloaking efficiency of ellipsoidal (a ^ b c) cloaks. 



shape 


direction 


inner 








outer 








rj 


ellipsoid 


a 


ai = 


0.1m, bi 


= 0.2 m, ci 


= 0.4m 


a2 = 0.2 m, &2 


= 0.4 m, C2 


= 0.^ 


^m 


99.55% 


ellipsoid 


b 


ai = 


0.1m, bi 


= 0.2 m, ci 


= 0.4m 


a2 = 0.2 m, 62 


= 0.4 m, C2 


= 0.^ 


^m 


99.56% 


ellipsoid 


c 


ai = 


0.1m, bi 


= 0.2 m, ci 


= 0.4m 


a2 = 0.2 m, 62 


= 0.4 m, C2 


= 0.^ 


^m 


99.52% 
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